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Abstract 

The low energy scattering of pions is investigated in the presence of electromagnetic 
interactions at leading order and at next-to-leading order for the amplitudes involving at most 
one pair of charged pions. The size of the electromagnetic corrections to the S-wave scattering 
lengths is found to be comparable to the size of the two loop strong interaction corrections. 
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1 Introduction 



Over the last few years, low energy tt — tt scattering has been increasingly recognized as provid- 
ing a perhaps unique window through which one might contemplate the actual mechanism of 
spontaneous breakdown of chiral symmetry in QCD (for recent discussions, see [|[]). It is nowa- 
days widely expected that this spontaneous breakdown is triggered by the formation of a large 
bilinear quark condensate < | gg | >~ —(250 MeV)^, which provides the main contribution 
to the masses of the light pseudoscalars mesons 

_ (rriu + mrf) < I gg I > 

However, on the basis of our present theoretical understanding of non-perturbative phenomena 
in asymptotically free gauge theories, the situation where the condensation of quark- ant iquark 
pairs in the vacuum is much weaker, say < | gg | >~ —(100 MeV)^, cannot be excluded 
(see for instance ||^). It turns out that a significant deviation of the ratio ( |1 . ID from unity 
influences low energy n — tt scattering 0] : The weaker the condensate, the stronger do pions 
interact at low energies. Unfortunately, the experimental informations available at present do 
not reach a precision which would allow to rule out any of the possibilities mentioned above. 
In this context, planed experiments such as the measurement of the lifetime of 7r+ — atoms 
at CERN 1^, or new high statistics Ki^ experiments by the E865 and KLOE collaborations, 
respectively at BNL and Da$ne 0, are of particular interest. 

The expected precision of these forthcoming data has triggered quite some activity at the 
theoretical level. In particular, the vr — vr scattering amplitude has been obtained at next-to- 
next-to-leading order both in the framework of generalized chiral perturbation theory |^ , where 
the ratio ( |1 . 1| ) is kept as a free parameter, and for the special case of standard chiral perturbation 
theory [^, where this ratio is exactly equal to one at leading order. In the former case, the 
whole range allowed by the present experimental value of the / = S-wave scattering length, 
Oq = 0.26 ± 0.05, is covered as the ratio (|1 . 1|) is varied between and 1. In the latter case 



the prediction, as taken from the recent literature P, 0, is in the range aj] = 0.205 — 0.217, 



depending mainly on the evaluation of the 0{p'^) counterterms (this question is presently under 
further investigation |jll[). These last numbers, when compared to the corresponding leading 
order |jT2[ and next-to-leading order values |T3|, Oq = 0.16 and = 0.20, respectively. 



strongly suggest that contributions beyond two loops are negligible (this conclusion also holds 
for the generalized case [0]). 

In order to assess the prediction of the standard far as low energy tt — vr scattering 

is concerned, and thus be able to detect any deviation from the strong condensate scenario 
( |1.1D once more precise data become available, it becomes mandatory to watch out for other 
contributions, which might compete with the two loop strong corrections discussed above. 
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For instance, the aforementioned results were obtained upon neglecting all isospin breaking 
effects, i.e. upon ignoring the quark mass difference m„ — Trid, and without taking into account 
electromagnetic interactions. In the present paper, we are interested in the evaluation of the 
corrections induced by the latter, which are expected to be the dominant ones (just as they 
dominate, for instance, the difference between the charged and neutral pion masses [p!5| , |1^], 
whereas 7^ 771^ only affects both the pion masses and scattering amplitudes starting from 
next-to-leading order). We shall work within the framework of chiral perturbation theory 



|T6| , |T^, 0, extended to include electromagnetic interactions [|T^, The central object of 
our analysis is the generating functional 



e 



^^(^^'^^'^'P'^^^'Q^^ = < 0; out I 0; in >.„a„.,,,Q,,Q, , (1-2) 



where the expectation value is obtained by integration over gluon, quark and photon config- 
urations, weighted by the action corresponding to the QCD Lagrangian with electromagnetic 
interactions and in the presence of the external sources v^, a^, s,p, Ql, Qr, 



^ = ^QCD + + (lY[v^,{x) + l5a^{x)]q - q[s{x) - i-f5p{x)]q 

+A, qY { Ql{x) {^^] + Qr{x) CA^] } q • (1.3) 



2 / ' \ 2 

In the above expression, C^qp^ stands for the QCD Lagrangian with massless light quarks, 
while £° stands for the Maxwell Lagrangian of the photon field A^. In both cases, some gauge 
fixing and the corresponding Faddeev-Popov ghosts are understood. The last term provides 
the interaction of the A^/ light quarks q = u,d,... with the electromagnetic field once the 
sources Ql{x) and Qr{x) are put equal to the constant charge matrix Q appropriate for Nf 
light flavours. The structure of the generating functional Z above has been obtained to one 
loop order in |2^ for Nf = 3 within a systematic framework which combines the usual chiral 
expansion [|16], 17, 14, 18| with the expansion in powers of the electromagnetic coupling e (see 



also [0). For our purpose, we need only to consider the case Nf = 2. However, for those 
aspects which are independent of the actual number of light flavours involved, we shall quote 
the results for arbitrary Nf. Electromagnetic corrections in the two flavour case within this 
systematic framework have also been considered recently in 123], and we shall compare 



with our present work in due course. Let us mention, however, that the only application that 



was discussed in [|^, ^ is the scattering of neutral pions, where photon exchanges between 



pions do not contribute before the two loop order. In the meantime, the formalism of ^ 
has been used to compute the corrections to the lowest order formula |^ for the lifetime of 
pionium in Ref. |^ (for other recent approaches to this problem, see [^). This computation 
obviously involves the scattering amplitude tt+tt^ —>■ tt^tt^, but, for the purposes of the bound 
state problem, at an off-shell point, whereas we shall consider the same amplitude for free, on- 
shell, pions. Whenever comparison was possible, we found agreement between our expressions 
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and the results of |25|. Finally, electromagnetic corrections to low energy tt — tt scattering 
within different contexts and/or frameworks have been considered previously in Refs. [P^] , 
p8| and [^]. Unless otherwise stated, we shall work within the framework of standard chiral 
perturbation theory, where is assumed to hold. 

The rest of the paper is organized as follows. In Section 2 we investigate lowest order 
electromagnetic effects on the vr — vr scattering amplitudes and discuss the S-wave scattering 
lengths. Section 3 is devoted to the construction of the generating functional to one loop 
order. We give a complete list of counterterms for Nf = 2 and for constant sources Ql{x) = 
Qr{x) = Qi and compute the corresponding divergent part of the generating functional. The 
amplitudes for ti^-k^ —>■ tt^tt^ and tt+tt^ tt^tt^ are worked out to one loop order in Section 
4. Scattering lengths are discussed in Section 5. In particular, the definition of the scattering 
length for tt+tt^ tt^tt^ requires a proper treatment of the infrared singularities which affect 
the amplitude. A final Section summarizes our results, and details on various technical aspects 
have been gathered in an Appendix. 



2 Lowest Order Electromagnetic Corrections 



It has been shown in |20] that in presence of electromagnetism, a consistent expansion scheme 



is obtained if the electric charge e is considered as a quantity of order p in the chiral counting, 

e, Ql, Qr ~ 0{p) , (2.1) 

where p denotes a typical momentum, much smaller than the scale Ah ~ 1 GeV at which the 
(non-Goldstone) hadronic bound states appear. At leading order in this combined expansion, 
and for an arbitrary number Nf of light quark flavours, the generating functional is given by 
the tree graphs constructed from the effective lagrangian [|19| 



-^Ff-'F^, - ^{d-Af + C{ QrUQlU^ ) . (2.2) 



For the notation, we follow Ref. |2^, where the transformation properties of the various 
quantities which enter ( p.2| ) under the chiral SU(A^/) xSU(A^/) symmetry group are displayed. 
Besides the scalar and the pseudoscalar sources, the field x contains also the mass term for the 
pseudoscalar mesons at leading order, x = 2B A4 H — ■, where M. is the (diagonal) quark mass 
matrix. The low energy constant B describes the bilinear light quark condensates, 

<0|mm|0> <0\dd\0 > 
B = = = ■■• (2.3) 
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whereas F is given by a two point correlation function of the vector and axial currents at 
vanishing momentum transfer, 

FS'^' = jj rf'x<0|T{(gL7'^A»(a;)(gR,7;.AV)(0)}|0> . (2.4) 

In both cases, | > stands for the vacuum in the chiral limit and in the absence of electro- 
magnetism, I >= I Vt > |_A4=o,e=o- The matrices A°/2 are the generators of the corresponding 
SU(A^/) flavour group. The covariant derivative ci^, defined as 

d^U = d^U - i{v^ + QrA^ + a^)U + iU {v^ + QlA^ - a^) , (2.5) 

and the last term of (|2.2|) contain the "spurious" Q1{x) and Q\{x), which play the role of 
sources for insertions, into the QCD Green's functions, of the electromagnetic vertex opera- 
tors A^qi^ (A'^/2)7'^gL and A^gR (A"/2)7'^gR, respectively. The low energy constant C gives an 
electromagnetic contribution to the charged pseudoscalar masses, for instance. 



Ml, = (m„ + md)B , 

Ml± =(m„ + mrf)5 + 2C-— , (2.6) 



which yields 



F2 

C_ ^ Ml^ - Ml, 
pi 2e2F2 



(2.7) 



For F = 92.4 MeVQ, this gives Z=0.8. Alternatively, C is given by the same correlator as in 
(p. 41), convoluted with the free photon propagator function, 

CS""' = ^ J d^xD^^x) <0|r{(gL7MA"gL)(x)(gR7.AV)(0)}|0> . (2.8) 

Finally, the penultimate term of (|2.2| ) acts as a gauge fixing. Notice that C is independent of 
the gauge parameter a, since the correlator in (|2.8|) , involving conserved currents, is transverse. 
Since this correlator itself is an order parameter of spontaneous breakdown of chiral symmetry, 
it has a smooth behaviour at short distances. This in turn allows to convert the representations 



41) and (p.8|) into sum rules |31, B3] via unsubtracted dispersion relations. Upon saturating 



these sum rules with resonances, this yields an independent evaluation of C which is compatible 



with the number given above [p2| , p!9| , p^ . For our purposes, the leading order determination 
( ^.7[ ) will be sufficient. 

^ Since in this Section we work at lowest order, we identify F with the pion decay constant = 92.4 MeV 
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In fact, for the two flavour case, the sphtting of the pion masses is the only direct effect 
induced by the presence of C. This is best seen in the so called a-model parametrization of the 
field f/, 



U = Jl-mi2F^ + tl with 0=1 1 , (2.9) 



F \ V^TT , — TT 



in which the last term of ( p.2| ) gives no interaction vertices for Ql{x) = Qr{x) = Q, with 
(5=exdiag(|, — |) the two flavour charge matrix, 

C{QUQU+) = 2C ■ 1^ ■ TT+TT- . (2.10) 

This means that if we introduce the "isotriplet" states | 7i"'{Ma,'p) >, a = 1,2,3 (we use the 
Condon-Shortley phase convention), 

\7t\M^±,P)> = -l=(\n+{M^±,p)> +\7i-{M^±,p)> ) , 

\7r\M^±,p)> = -^(|7r+(M,±,p) > - |7r-(M,±,p) > ) , 

\tt%M^o,p}> = |7r°(M^o,p) > , (2.11) 

the amplitudes A"'^'''^'^{s,t,u) for the processes tt^tt^ tt'^tt'^ are expressed, after subtraction of 
the one photon exchange Born terms, by means of a single amplitude, 

A'''''"^is,t,u) = { A{s\t,u)S''^5"^ + perm} + ©(eV, e^) , (2.12) 

where, as in the absence of electromagnetism, 

. / , N s — 2mB 

A{s\t,u) = (2.13) 

at lowest order, with m = |(mu + m^). Several remarks are in order at this point. First, 
the observations contained in Eqs. ( |2.1CI| ), ( |2.12| ) and ( p.l3| ) were made previously by several 
authors (see e.g. |34| and references therein). Second, we have added a contribution 0(e^p^, e^) 
in ( p.l2| ) to indicate that this representation might not be maintained once next-to-leading 
electromagnetic corrections are taken into account, even after suitable subtraction of the long 
range one photon exchange contributions. We shall come back to this point in the following 
Section. Finally, at first sight the representation ( p.l2|) suggests that the leading electromag- 



netic corrections induced by the last term in ( p.2| ) do not lead to isospin violations in tt 



vr 



scattering. However, and as explicitly indicated by the notation in (|2.11|) , the "isotriplet" states 



have unequal masses, and these mass differences lead to isospin breaking contributions to the 



5 



amplitudes (for instance, the thresholds of various reactions can be distinct). Considering the 
S-wave scattering lengths, we write 



ao(00;00) -- 
ao(+0;+0) 

«o(+-;00) 

ao(H — 
ao(++ 



3V"'0/'str -r 3 v"'0/'str + Aao(00;00) 
^{a^o)str + Aao(+0;+0) 

" + Aao(+-;00) 
f Aao(+-;+-) 



3(0-0)8^ 



3(0-0)8^ 



0/str 



++) = {al)str + Aao(++;++) 



(2.14) 



where (ao)8tr and (ao)str denote the strong isospin / = and 1 = 2 S-wave scattering lengths in 
the absence of electromagnetic corrections, and Aao{ab;cd) are the corresponding corrections 
induced by electromagnetic effects. At lowest order, one has (ao)8tr = 7M^/327rF^ and (cio)str = 
— M^/lGvrF^, where M^. denotes the leading order pion mass for e = 0, i.e. = 2mB. From 
the expressions ( |2.6D , M^^ should therefore be identified with M^o at this order. However, 
it has become customary to quote the values of the lowest order S-wave scattering lengths, 
obtained by Weinberg more than thirty years ago 0], with the value of the charged pion mass 
M^± = 139.57 MeV |5D| assigned to M^, i.e. 



(^o)str 



7M\ 
327rF2 



0.16 



167rF2 



-0.045 



(2.15) 



Adopting this definition, we obtain the following values for the shifts in the S-wave scattering 
lengths (A, = M2± - M^o) : 



Aao(00;00) = 

Aao(+0;+0) 

Aao(+-;00) 

Aao(+-;+-) 

Aao(++;++) 



327rF2 
A. 

327rF2 
A^ 

327rF2 
167rF2 
167rF2 



(-6.4%) 
(+6.4%) 
(-2.1%) 
(+6.4%) 
(+6.4%) . 



(2.16) 



The absolute variation is at most A^/167rF2 = 0.003. The relative variations are shown between 
parentheses. Furthermore, the scattering lengths of the four amplitudes which involve at least 
one pair of charged pions satisfy the usual isospin relations. In order to make this apparent, 
we introduce the modified scattering lengths 

_ A. 



2-0)811 



''O/str 



+ 



327rF2 

A^ 
167rF2 



0.166 



-0.042 



(2.17) 
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and obtain 



ao(00;00) -- 

ao(+0;+0) 

ao(+-;00) 

aoi+-;+- 

ao(++;++ 



1^0 



2^2 



2"0 

„Uq -t- „Uq 



A. 



87rF2 



1^0 



1^2 



an + 7?a 



(2.18) 



Thus, at leading order, we can distinguish two consequences of electromagnetic interactions for 
the n — 71 scattering lengths : A shift in the strong S-wave isospin scattering lengths (ao)str, 
described by Eq. (|2.17|) , and an additional explicit isospin breaking contribution to ao(00; 00) 
in Eq. ( [2.181 ). These corrections to the scattering lengths are comparable in magnitude to the 
pure strong interaction two loop effects as evaluated in the recent literature p, ^ . 

The expressions ( |2.16|) above disagree with those obtained by Maltman and Wolfe in Ref. 
p9| . Of course, on general grounds |^ the result ( p.l6|) does not depend on whether one uses 
the parametrization (p.9|) or the exponential parametrization 



U 



(2.19) 



as Maltman and Wolfe do. Rather, these authors have not taken into account the isospin 



violating contributions that come from the first term of ( p.2| ) via the mass difference between 
charged and neutral pions. Therefore, their expressions for the leading order electromagnetic 
corrections to the scattering lengths are not correct. 

Interestingly enough, at the same order, the relation ( |2.12| ) also holds in the generalized 
case0, and with the same expression ( p.l3| ) for A{s\t,u). The difference with the standard case 
comes from the fact that 2mB no longer represents, for e = 0, the only contribution to the pion 
mass even at leading order. Instead, one has 



2mB 



a 



(2.20) 



where the parameter a (which is not related to e^/Air !) varies between 1, the standard case 
of a strong condensate, and 4, the extreme case where the condensate would vanish^. Corre- 
spondingly, the leading order S-wave scattering lengths become 0], for e = 0, 



la, 



O/str 



967r F2 



(5a + 16) , 



(^o)str 



487r F2 



(a 



4) 



(2.21) 



^J. Stern, private communication. 

■^Within the SU(3) framework, a can be related to the quark mass ratio ms/rh : a = 1 corresponds to 
the standard value ms/fh = 25.9 at leading order |l^, whereas for a — A, the quark mass ratio drops to 
TOs/m = 6.3 H. 
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However, the result ( p.l^ ) for the scattering lengths holds as it stands, provided the definition 
of the corrected isospin scattering lengths is suitably modified to read 

a'o = iOstr + i (4 - «) , 

al ^ {alUr + I (4 - a) . (2.22) 

As the condensate becomes weaker {i.e. as a grows towards 4), the magnitude of the correction 
induced by electromagnetic effects in the scattering lengths Oq decreases, while the explicit 
isospin breaking component of ao(00; 00) is not affected by the size of the condensate. 

We conclude therefore that electromagnetic corrections to the S-wave scattering lengths 
are sizeable at leading order. Whereas they remain small as compared to the one loop strong 
interaction corrections, they are of the same order of magnitude than the two loop corrections 
0, H, |1^. In the sequel, we shall investigate how these numbers are affected by next-to-leading 
corrections. Our study will henceforth be limited to the standard case. 



3 The Generating Functional to One Loop 

In this Section, we construct to one loop order the generating functional Z{v^, a^, s,p, Ql, Qr) 
in the presence of electromagnetic interactions. 

At next-to-leading order, this generating functional involves one loop graphs with vertices 
from and tree graphs with vertices from C^'^^ and at most one vertex from the next-to- 
leading effective lagrangian'^ 

The first term contains the purely QCD low energy interactions among the pseudoscalar mesons. 
In the case of two fiavours and for e = 0, its expression was given in Ref. [|14|. In the presence 
of electromagnetic interactions, it readsQ : 



■*This decomposition of £^'*-' corresponds to the minimal number of sources Ql,r involved. Terms proportional 
to the electric charge e are also present in Cpi via, for instance, the covariant derivative d^. 

^We use here the SU(2)xSU(2) formalism, rather than the 0(4) formalism of jlj], particular to the two 
flavour case. We also omit the Wess-Zumino term, which describes anomalous couplings between photons and 
pions (for a review, see ]3^). 
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+ h ( G^,UG^'"'U+ ) + y ( Gf^.d^Ud^U^ + G^J^'U+d'^U ) 



+ \{h,- h) Reidetx) - h ( G'l^.G^'^^ + G^^G'^''^ ) . (3.2) 

For the definition of tlie quantities G^i,^, we refer the reader to Eqs. (A. 6) and (A. 13) of the 
Appendix. 

The loop graphs with vertices from C^"^^ generate divergences, which are absorbed into 
the renormahzation of the low energy constants of C^^\ In order to compute these divergent 
pieces, one may evaluate, for an arbitrary number of space-time dimensions d, the functional 
determinant of the quadratic part of the quantum fluctuations of the pion and photon fields 
around their classical configurations in the presence of the external sources (see Ref. pO| and 
the Appendix). For our purposes, we need only to consider the two flavour case (the relevant 
expressions for Nf arbitrary can be found in the Appendix), and with the sources Ql{x) and 
Qr{x) put to their constant value, 

Ql{x) = Qr{x) = Q , Q = ex diag(|-i) . (3.3) 

Furthermore, unless otherwise stated, from now on we restrict ourselves to the Feynman gauge 
a = 1. Upon using the identity 

(Q)' = liQ') , (3.4) 

we then obtain, from the expressions (A. 9) and (A. 14) of the Appendix, 

1 1 



.'one loop, div 



167r2 d-A 



d X 



I ^ ( d^U^d^U f + ^ {d^U+d'-U ) ( d^U+dM ) 

-- ( GtUG^^'''U+ )--{ G^d^Ud^'U^ + GLd^U+d^U ) 
6 ^ 6 ^ ^ 

+l(x^x)-^{ G^.G^"^ + G^^G'^^'^ ) + Reidetx) 



30 
27 Z 

20 
3F2 



+ -) F2 ( df'U+d^U ) ( ) + 2ZF^ ( d''U+d,,U ) ( QUQU+ ) 
( d^f/+Qf/ ) ( d^U+QU ) + ( d'^U+QU^ ) ( d^UQU^ ) ) 



4 

+2ZF2 ( d^'U+QU ) { d^UQU+ ) 
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+ (] + 2Z] { x^U + f/+x ) ( QUQU+ ) 



4 

+ Q - ( (xf/+ - Ux^)QUQU+ + ix^U - U+x)QU+QU) 



+^ ( d.U+MQ), Q]U + d^U[{clQ),Q]U+ ) 

3 4,„ „ +,2 



+ [- + 3Z + UZ'j F" ( QUQU+ ) 

The covariant derivatives {c^iQ) are defined in the Appendix, Eq. (A. 12). 

The complete list of possible counterterms in £e2p2 and was given for A^/ = 3 in Refs. 
pO| , PI) . For Nf = 2, the number of possibilities decreases, due to additional trace identities 
for products of 2x2 matrices. In the simpler situation described by (|3.3|) and (^.4|), we find : 



and 



= {k,{d^U+d^U){Q') 

+k2{d''U~^d,,U){QUQU+) 

+h ( ( d^U+QU ) ( d^U+QU ) + ( d^UQU+ ) ( d^UQU+ ) ) 
+h { d^'U+QU ) { d^UQU+ ) 

+h{x^u + u+x){Q') 
+h{x^u + u+x){QUQU+) 

+kr { (xf/+ + Ux^)Q + ix^u + u+x)Q ){Q) 

+ks { (xf/+ - Ux^)QUQU+ + (x+f/ - U+x)QU+QU) 
+kg { d^U+[{c^j,Q), Q]U + d^U[{clQ), Q]U+ ) 
+fcio((4Q)f/(cLMg)f/+) 

+ku ( (cnQ) ■ (crQ) + {clQ) ■ (clQ) ) } , (3.6) 

C,, = F^ { ki2 ( )' + fcis ( QUQU+ ){Q^) +ku{ QUQU+f } . (3.7) 



In the large Nc limit, the constants ki are of order 0(1) for i = 1, ...11, and of order 0{1/Nc) 
for i = 12, 13, 14. One combination of these constants is Zweig rule suppressed, 

k5 + ke- ~ 0{1/Nc) . (3.8) 
5 
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The renormalization of the low energy constants of is given as 



hi 



1,... 7, 
1, 2, 3 , 
1,... 14 



with 



,d-i 



X 



^ ( -ihi47r + r'(l) + ll 

167r2 \ d-4: 2^ ^ ' ^ 



The coefficients 7j and 6i were computed in [|IJ], and can be read off from Eq. 
addition yields 



(3.9) 

(3.10) 
]5D, which in 



(Ti = 
(74 = 

ar = 

0"lO 
0"12 



27 _ 1 7 
20 5 ' 



2Z, 
0, 
= 0, 

_ 3 _ 
2 



(T2 = 
0-5 = 

as = 

(^11 

84 72 ^ 
25^ ' ^13 



2Z, 

4 5"^' 



0, 



-3 - |Z 

5 



12 72 



2Z, 



1 _ 

4 

1 

4' 



I + 3Z + 12^2 . 



(3.11) 



Up to a reshuffling of the indices^ of the low energy constants ki, the list (|3^ ) and (|3.7| ) 



coincides with the one given by the authors of Ref . |2^ , except for the counterterm multiplied by 
kj, which they have omitted. This term, however, cannot be rewritten as a linear combination 
of the remaining structures that appear in Ce2p2, and must therefore be included, even if fcy 
turns out not to be renormalized. In order to illustrate this point, we consider isospin breaking 
in the light quark condensates, which appears only at next-to-leading order, but receives no 
contributions from the pion loops. 



< Q\uu — dd\fl > = ^{rrict 



(3.12) 



For e = 0, the result involves only the "high energy" constant /13 [Tj], which is a convention 
dependent short distance subtraction. In the chiral limit (or, to that order, for rriu = irtd) and 
with e 7^ 0, the isospin violation in the condensates is given solely by the constant kj. Both /13 
and kj are finite, 6^ = a-j = 0, but, in contrast to the former, kj is an order parameter for the 
spontaneous breakdown of the SU(2)xSU(2) chiral symmetry. 

As far as the remaining /cj's are concerned, Ref. also shows some discrepancies with 



our values of the renormalization factors cij. In a subsequent version ||2^, its authors brought 
their findings into agreement with our results ( |3.11| ). Finally, the divergent part ( |3.5| ) contains a 



renormalization of the Maxwell Lagrangian of the photon, which was not mentioned in E^, 



°We have numbered our kiS, such as to make the correspondance with the SU(3) counterterms Ki oi \ 2Q , 2 
easier. 
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As a second application, we also give the expressions of the pion masses at next-to-leading 
order, 



Ml. = 2fflB{l + H||«A)+e'r.,.W + ^l.>(^)-3||5ln 

8 

-2—{md - niuf h - ^B{md - m„)e^ kj 
+2e^F^ ( ^(1 + ^) + e^/C:4/i) - (3 + 4Z)-^ In 



g 

■-B{md-mu)e^ k7 (3.13) 



with 



20, 
20, 



/C;± = --[kl + k;-kl-l{23kl + kr + 18kl) 



/C:± = -H[2Z(fc[ + fc^')-^^[3-fcu] • (3.14) 



As expected, these expressions turn out to be independent of the choice of the renormalization 
scale fi, which provides a non trivial check of the relations ( p.ll| ). The low energy constants 



I7 and ^7 induce corrections of the type {rrid — and (m^ — mu)e^, respectively. While 

the latter do not contribute to the difference M'^± — M^p, the former have been estimated in 



15, 14] and were found to be negligible as compared to the leading order electromagnetic mass 
difference (p^ ). If, on the basis of naive dimensional analysis, we assume an upper bound on 
the various constants fc[(/i) at the scale /i = Mp, 

I fc[(M,) I < ^ , (3.15) 

we obtain 

^\lCAMp)\ < 0.9x10-2 

^\lCl±{Mp)\ < 1.5x10-2 

e2|/C:±(M,)| < 0.3 X 10-2 , (3.16) 

so that the corrections of order 0{e^m) and 0{e^) to the pion masses turn out to be numerically 
small. They are, for instance, comparable to the error on the 0{m?) corrections generated by 
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the uncertainty assigned to /s in 0, viz. I3 = 2.9±2.4 or {2M^o/F'^) ll{Mp) = (0.4±1.6) x IQ-^. 
More elaborate statements require more information on the values of the constants ki. In the 
present paper, since we are mainly interested in assessing the expected order of magnitude of 
electromagnetic corrections, we shall be content with the cruder estimate (|3.15|) . 

Before we turn to vr — vr scattering, we should like to stress that the result ( |3.11| ) only holds 
in the Feynman gauge a = 1. In general, the low energy constants ki, and thus their divergences, 
depend on the choice of gauge. This gauge dependence has been displayed explicitly for some 
of the SU(3) counterterms Ki of |^ |21| in Ref. |3^, but a systematic investigation has not 
been undertaken so far. Such a calculation would, for instance, be useful to check that the 
expressions for physical quantities, like the pion masses above, are not only scale independent, 
but indeed also gauge invariant. 



4 The One Loop n — n Scattering Amplitudes 

In this Section, we compute the amplitudes for tt^tt^ tt^-k^ and for tt+tt^ 71^71^ scattering 
at next-to-leading order in the presence of electromagnetic interactions and for niu = rrid. The 
amplitudes involving charged pions only are however not considered here. 

For the process vr'^vr'^ tt'^tt^, the only differences with the calculation in the absence 



of electromagnetic interactions |]14[ consist in taking into account the contributions from the 
counterterms of Ce^p^ (those of do not contribute here) and in keeping track of the masses 
of the pions in the internal lines of the loop graphs. The amplitude is then given by 

+ 48^?^ (^1 + 2/^2 - 3L. - 3) [ - 2M2o)2 + (t - 2M'^or + (u - 2Mlpf ] 

[{s - Ml,fj^_{s) + {t- Ml,fj^_{t) + {u- Mlpfj^_{u) ] 
+ ^,[Jm{s) + Jm{t) + JM] , (4.1) 

/coo ^ (3 ^ _ J^^^ _ 3^^ _ ^zk, , (4.2) 



where we have set 



with 



and 



32^^ + 



= I k^ + ln-f] , z^7, (4.3) 



In T# • (4-4) 
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The loop functions Jpq are as defined in (see also the Appendix). The subscript identifies 
the charges, and therefore the masses, of the two pions in the loop. As observed in [^, the 
contribution of the electromagnetic counterterms can be eliminated upon replacing F by the 
neutral pion decay constant -F^o[], 



which describes the coupling of 7r° to the axial current, 

< ^ I ^(w7m75M - (^7m75'^)(0) 1 7r°(p) > = iy^ F^o , (4.6) 



for TUu = nid- However, since there is no accurate determination of this matrix element 



this substitution is of little avail in practice. Instead, we shall use the charged pion decay 
constant F^^ in the absence of electromagnetism, defined as 

< n I (M7^75rf)(0) 1 7r~(p) > |e=o = i V2p^ F^ , (4.7) 

with the value Fj, = 92.4 MeV |39[ as extracted from the 7r£2(7) decay rateQ, and the 

relation to F given by [T^ 



In the case of the process tt^{pi)tt^{p2) —>■ 7r°(A;i)7r°(A;2), one has to consider in addition 
loops with virtual photons. These affect the wave function renormalization of the charged 
pions and introduce a long range component into the scattering amplitude through the vertex 
correction graph of Fig. 1. With the Condon-Shortley phase convention, the amplitude reads 

o _ /If 2 

A+-''\s,t,u) = - — ^ - B{s,t,u) - C{s,t,u) , (4.9) 




Figure 1 : The one photon exchange electromagnetic correction to the strong vertex. 



^The coefBcient of Lti differs by a factor of two from the corresponding expression in |22 . 
^This extraction itself is not without ambiguities. In order to fix these, a more extended framework, aUowing 
for a full treatment of radiative corrections in the semileptonic decays of the pseudoscalar mesons, is required. 
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where B{s,t,u) represents the unitarity corrections 



B{s, t, u) 



F4 



-e'F' (2s - 4M^±) G'+_^(s) | 



t - 2M^± + -f 

V 



12F4 



^2 



u - 2M^± + ^ 
u 



(4.10) 



whereas the contributions from O(p^) tree and tadpole graphs have been collected in C{s,t,u) 
(E, = M^^ + M^o), 

C{s, t, u) = 1^ I - ^ - 4A. - ^ [AMt. - IMl^Ml. + ^M^o 



2 + (3 + —)ki + -^k2 + 3h + ^Zki - 6L^ + 41n — ^ 



9 



9 



M% (ML 10M2 L 



327r2F4 



9 



( 2M,2± - M^^o ) + M^oh 



8Zk2 + Sks + AZk^ - 2(1 + SZ)kQ - (1 - SZ)k^ 



967r2F4 



1 



(1 + ^) {^^K-2Ml,Ml^L^) 



--M^o)s' + M^± (i" + u') 



487r2F4 



(t - S,)2 + (m - E,)^ 



(4.11) 



The loop function G+_^(s) involves three propagators, and comes from the Feynman 
graph of Fig. 1, 



G'+_-y(s) 



— z 



1 



(27r)4 {q^ -m'l){q^ -2q-pi){q'^ + 2q-p2) ' 



(4.12) 



While being ultraviolet finite, this function develops an infrared singularity as the photon mass 
is sent to zero. In order to display this singularity, we express the function G+_^(s) in 
terms of the dilogarithm or Spence function Li2(x) (see the Appendix). For s > AM^±, and for 



15 



m-y ^ 0, it is given as follows : 



is) = 



where 



327r2 
+2 



sa 



In 



-In I 



Ati^ ^ ^ fl-a 

h InM 

3 \l + a 



+ 21n(cT) 



In ( I ] + m 

1 + (7 



'1- 



4M2, 



}, (4.13) 



(4.14) 



Both the real and imaginary parts of G+_-y(s) diverge at threshold. Another infrared singular 
piece, coming from the wave function renormalization of the charged pions, appears in C(s, t, u). 
As is well known, in order to obtain an infrared finite cross section, one has to consider the 
processes with emissions of soft photons. This question will be addressed in the next Section. 
For the moment, we furthermore notice that both B[s, t, u) and C(s, t, u) contain terms which 
involve inverse powers of the Mandelstam variables t and and whose coefficients vanish in 
the limit where the neutral and charged pion masses become equal. In fact, the coefficient of 
these terms in C(s, t, u) is rather small : Expanding the logarithm in powers of A^, one has 



2M\Ml±L^ 



m;Jo / A^ 



1 + O 



( A. 



(4.15) 



which yields a negligible contribution even near threshold, where t u —^n- Keeping only 
terms which are at most of order A^/M^o, the expression of C(s,t,M) simplifies somewhat 



C(s, t, u) 



Ml. 



327r2F4 
1 



1 

487r2F4 



-4M,2o - 7A. - m + e^F W° 



4e^F2ln. 



mi 



29 ^ 
18 Ml, 
2 A! 



Ml. (/- - 1) - 3 A. - ^ + e^F^/Cf 



3 4M2, 3M;J„ 
5 1 A. 1 A! 

I'2 — - — - TTT- + - 



6 2 Mlo 6 Mto 



{s-2Mlo){s-2Ml^) 



{t - E,)2 + {u- E^f 



(4.16) 



with 



±0 



/C^° = 8ZA;2 + 3A;3 + 4ZA;4-2(1 + 8Z)A;6- (1 -8Z)A;8 . 



4^ - 32^- 
(3 + —)ki + -—k2 + 3A:3 + 4Zh - 
9 9 



(4.17) 



In the sequel, we shall always neglect terms which are smaller than A^/M^o or than Ajr/M^o- 
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As far as the representation ( 2.12| ) is concerned, if it were still correct at next-to- leading 
order, it would entail the following relation among the two amplitudes we have obtained above : 



- = + + , (4.18) 

once the one photon exchange and infrared singular contributions have been subtracted from the 
amplitude A'^~''^^{s, t, u). Even without figuring out how to perform this step in an unambiguous 
way at the present stage, it is clear that such a relation cannot hold : Not only are the unitarity 
corrections due to two pion intermediate states different from what the relation ( |4.18| ) would 



require, but also the contributions of the electromagnetic counterterms at order 0{e m) do 
not satisfy this relation. If one considers in addition the scattering amplitudes involving four 
charged pions, one can easily see, even without doing a full calculation, that they will involve 
the constant ki4^ of Cgi, which does not contribute to the amplitudes involving neutral pions 
or to the masses. Therefore, we conclude that the representation (|2.12| ) in terms of a single 
amplitude A(s|t, u) is violated by next-to-leading corrections both of order 0{e^p^) and of order 
O(e^). The meaning to be assigned to the amplitude defined by Eqs. (22), (23) and (24) of 
Ref. 1^ thus remains unclear to us. 

More than twenty years ago, Roig and Swift |^ had already investigated electromagnetic 
corrections of order O(e^p^) to the vr — tt scattering amplitudes. Their approach was not 
systematic from the point of view adopted in the present article : the pion loops and the pion 
mass difference were not taken into account. Thus, up to the Born terms, only the photon loop 
induced wave function renormalization of charged pions and vertex correction graphs of the 
type shown in Fig. 1 were included in their calculation. Their claim, however, was that the 
resulting amplitudes were all finite. We have checked that this is not the case for tt+tt^ tt^tt" 
(the same conclusion was reached independently in Ref. [^]). Working within the framework 
adopted by Roig and Swift, we obtain (now in an arbitrary gauge a) 



Ats''\s,t,u) = -— ^|l + Y^-2e^(.-M^^,)G+_,(.)-e^J+_(^) 

-Ge^Aj. (4.19) 



167r2 



31n^ _21n^ 



Thus, although we end up with a gauge invariant result, it however diverges as d 4. This 
disease is cured by the counterterms fcj. We have also checked that the amplitudes involving 
four charged pions are divergent in the framework of Ref. . Thus, omitting the contributions 
of the electromagnetic counterterms to the scattering amplitudes does not lead to a consistent 
result. 
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5 Scattering Lengths 



Evaluating the amplitude ( [4.1| ) at threshold, s = 4M^o, t 
expression for the scattering length, 



u 



0, we obtain the following 



ao(00; 00) 



327rF2 



1 + 



327r2F2 



8/1 + I6/2-3/3 + 4/4 + I3 + 18— f + 



A2 



+ 



327r2 



00 



(5.1) 



where we have neglected terms of order A^/M^o and beyond. At next-to- leading order, 
coincides neither with M^± nor with M^o for e 7^ 0. We shall define the correction Aao(00; 00) 
as in Section 1, Eq. ( |2.14| ), but now with 



7M\ 5 



327rF2 



1 + 



167rF2 



1 



847r2 F2 
1 M2 



- - 3- 21- 21 

/i + 2/2 /s H h H 

8 10 8 



127r2 F2 



fi + 2r2 - Ih 



3- 3 

-^4 + - 

2 8 



0.20 ±0.01 (5.2) 
-0.043 ± 0.004. 



For the numerical evaluations, we use the central values of the low energy constants k that were 
Mji = -2.3 ± 3.7, I2 = 6.0 ± 1.3, [3 = 2.9 ± 2.4, h = 4.3 ± 0.9. With the estimate 



given m 



( ^ ), we find 



p2F2 

^ -'^TT 1^00 

M2o 
e^F2 



±0 



1.0 ±0.9 



1.8 ±0.9 



0.5 ±2.2 . 



and 



Aao(00;00) 



-3.2 ±0.1) X 10" 



(5.3) 



(5.4) 



which corresponds to a relative variation of —8% as compared to the one loop value ao(00; 00) |str 
= —0.039 obtained in the absence of electromagnetism. We see also that the contribution of 
the counterterms fc[ to the correction, which is reflected by the uncertainty in (|5.4| ), is very 
small, and would remain so even if the estimate ( ^.15| ) turned out to be wrong by a factor of 
two or three. 

As we have mentioned earlier, the amplitude for tt+tt^ tt^tt^ contains an infrared 
divergent contribution generated by the long range electromagnetic interactions between the 
incoming charged pions and by their wave function renormalization. The threshold expansion 
of the real part of the above amplitude takes the following form {q denotes the momentum of 
the charged pions in the center of mass frame) 



ReA- 



-;00| 



S, t, u) 



AM\-Ml, e2 M,± 



F2 



e 
16 



+ ReA 



;00 



thr 



+ 0{q) 



(5.5) 
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with 



Re A 



thr 



-;00 



327r 



'-(CtoJstr + -laoJstr 



A p2 i\/r2 



A. 



487r2Fi 
A! 



^ [M2o(1 + 4h + 3h - 12U) - QFy {10 - Kf) 
[212 -40^1- 15/3 + 180[4] . 



For the amphtude of the reaction tt^tt^ —* tt^tt^, the corresponding expressions read 



ReA 



+0;+0, 



S, t, u) 



ReA 



+0;+0 
thr 



+ 0{q) 



(5.6) 



(5.7) 



and 



ReA 



+0;+0 
thr 



167r (a^)str + ^ 



A! 



327r2F2 

3^3 



2 + /C^'' + /C2 



12L 



2887r2Fi 



[9 + 4/1 + 8/2 

[ 10 - 24^1 - 48/2 + 9/3 + 36/4 ; 



(5. 



Numerically, we find 
1 



327r 
1 



ReA 



;00 



thr 



Ojstr 



Oylstr 



-1.2 + 0.7) X 10" 



ReA 



+0;+0 
thr 



- - a, 



0/str 



;i.3 + 0.4) X 10" 



(5.9) 



327r 2 

The error bars refiect the naive dimensional estimates (|3.15| ) of the contributions from the 
low energy constants /c[(Mp), which again are rather small. Although the infrared singularity, 
which contributes only to the terms of order 0{q^) or higher, affects neither the constant terms 
i?ev4+;'°° and i?ev4+r^° 

nor the long range part of the one photon exchange graph, scattering 
lengths ought to be defined from infrared finite observables. In the remaining part of this 
Section, we shall perform this last step in the case of the process tt+tt^ — > tt^tt^. 

An infrared finite cross section is obtained if one takes into account the emission of soft 
photons with energies below the detector resolution. At the order we are working, it is enough 
to consider the emission of a single photon with the O(p^) 7r"'"7r^7r'^7r'^ vertex taken from 
The corresponding amplitude reads 



A-" 



;007 



F2 



{ki + k2 



ML 



(2pi - k 



(2p2 - k 



7/^ 



ml - 2pi ■ k^ mi - 2p2 ■ k^ 



(5.10) 



where k-y and e^(/c^) are the photon momentum and polarization, respectively. The infrared 
finite cross section is then given by the sum 



atot(s;AE) = a+-'^\s) + a 



;007, 



s-AE) 



(5.11) 
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with 



cr 



+ -;00, 



A 



+-;00 



(5.12) 



being the cross section without emitted photon, whereas the cross section with one emitted 
photon reads 



a 



-;007| 



s;AE) 



a 



-;007| 



s;AE) 



GAn^s \s- AMl± 



1 + ^ In 



sa 



1 - a 



In 



(5.13) 



The infrared finite piece '°°'^(s; AE) involves an integration over the energy of the emitted 
photon, which, in the present case, can be done analytically, and reads 



cr 



-;007 



1 



647r3s F^ \ s- ml 



4:M\ 



1 + 

sa 



1-a 
l + a 



+ {s- M^o 



2 \2 



In 2 + — In 
2a 



a 



l + a 



[F{s;AE)~F{s;0)] 
(5.14) 



The function F{s; AE) is given in the Appendix. The maximal energy AE of the undetected 
photon is limited by the detector resolution AE^et or, for s close to threshold, by the available 
phase space (in the present case, the maximal photon energy at threshold is A^r/M^i = 9 
MeV), 



AE 



, s - AM', ^ ^ 
min ( . , AEdet 



2v^ 



(5.15) 



As seen from the above expressions, the infrared divergent pieces cancel in the sum ( |5.11| ). 

Following Roig and Swift ||2^, we may define the scattering length ao(H — ; 00) from the 
threshold expansion of the infrared finite cross section atot. 



O"tot 



s - AMI, 



327rs \s- AMl± 



AM\ - M\ M^± 



F3 



16 q 



+ 327rao(+-;00) + 0(g) . (5.16) 



Notice that we have not expanded the ratio of the phase space by the fiux factor in front of 
c"tot- Up to terms which are of higher order in the chiral expansion, the scattering length thus 
defined reads 



327rao(+-;00) = ReA^^^ 



+-;00 



In 2) 



-M\ - -A^ 
4 2 



(5.17) 



and does not depend on AE [^. However, as compared to Eq. (5.6), it receives an additional 
contribution, the second term on the right hand side of (5.17), which is specific to the infrared 
finite observable cxtot that we have considered here. Numerically, this represents only a tiny 
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^ 280 300 320 340 

Figure 2 : From bottom to top, the pure strong interaction cross sections for tt+tt^ —>■ tt^tt" 
at leading and at next-to-leading order, and the infrared finite radiative cross section a^ot (in 
millibarn) as a function of the center of mass energy (in MeV). 



modification, and we obtain 

Aao(+-;00) = (- 1.2 ± 0.7) x 10"^ , (5.18) 

corresponding to a relative variation between —2.3% and —0.6% as compared to the one loop 
value ao(H — ; 00) = —0.082 for e = 0. Upon comparing Eqs. (|5^ ) and ( ^.181) with the values 



obtained at leading order O(e^) in Section 2, Eq. ( |2.16| ), we conclude that the electromagnetic 
shift in the scattering length ao(00; 00) including the 0(e^p^, e"^) contributions is twice as large 
as the leading order O(e^) effect, whereas the variation of ao(H — ; 00) is not strongly affected 
by next-to-leading order electromagnetic corrections. 

Finally, we compare the total cross section crtot to the cross section for tt+tt^ tt^tt^ 
scattering in the absence of electromagnetic interactions. The cross sections are shown in Fig. 
2. We have taken a typical detector resolution of AE = 20 MeV. The thickness of the upper 
curve shows the uncertainty induced by the estimate ( |3.15 ) of the counterterms of CfApZ. The 



error coming from the uncertainties of the strong counterterms /i,2,3,4 are not shown. In Fig. 3, 
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we show the ratio of (Jtot and of the one loop strong cross section with the leading order cross 
section for e = 0. As one approaches the threshold, the long range Coulomb contribution and 
the flux factor make the cross section diverge. Below 300 MeV, these two effects become more 
important than the strong one loop corrections themselves. 




280 



300 



320 



340 



Figure 3 : The relative corrections to the leading order cross section with only strong inter- 
action corrections (lower curve), and with electromagnetic effects included at next-to-leading 
order (upper curve), as a function of the center of mass energy (in MeV). 



6 Summary 

In the present paper, electromagnetic corrections to low energy tt — vr scattering were 
investigated at leading and next-to-leading order, within a systematic and consistent framework 



20, 21] which includes loops with virtual photons and combines the usual chiral expansion with 



the expansion in powers of the electric charge. 

At leading order, the only direct effect of the electromagnetic interaction between quarks 
is the splitting between the masses of the charged and neutral pions. In particular, a single 
amplitude allows to describe all the different n — tt scattering channels at this order. All 
the S-wave scattering lengths for the channels involving at least one pair of charged pions 
can be expressed, as in the absence of electromagnetism, in terms of two suitaby redefined 
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isospin scattering lengths Oq and Oq, Eq. ( p.l8|) . Only the expression for ao(00;00) receives 



an additional contribution, which explicitly breaks isospin symmetry. We have found that 
numerically these corrections are comparable in magnitude to those induced at two loop level 
by the strong interactions alone, as evaluated in Refs. . 



In order to obtain the vr — vr scattering amplitudes at next-to-leading order, we have com- 
puted the generating functional in the presence of virtual photons to one loop in the two flavour 
case. Its divergences (in the Feynman gauge) and the complete list of relevant counterterms 



were given exphcitly and compared to the results obtained recently in Refs. |23[ . Applica- 
tions to the isospin breaking in the light quark condensates and to the pion masses were briefly 
discussed. Upon using naive dimensional estimates for the electromagnetic counterterm, only 
small corrections to the situation for e = were obtained for the latter. 

At next-to-leading order, the complete expression of the scattering amplitude with a pair 
of charged pions has been worked out. The counterterms are required both to absorb the 
divergences produced by the virtual photon loops, and to provide a gauge invariant result. 
Actually, the second part of this statement was not checked explicitly. It would require that 
the gauge dependence of the various low energy constants were analysed in detail. The S-wave 
scattering lengths for the processes tc^ti^ tt^tt^ and tt+tt^ tt^tt" have been computed. 
For the latter, a proper definition requires that one considers the process with the emission 
of a single soft photon, in order to obtain an infrared finite cross section. The results for the 
electromagnetic shifts in the scattering lengths at next-to-leading order, 

Aao(00;00) = (-3.2 ± 0.1) x 10"^ 
Aao(+-; 00) = (-1.2 ± 0.7) x 10~^ , 

confirm the conclusion reached at the end of the leading order analysis : Although these 
numbers are small as compared to the corresponding values of the one loop strong scattering 
lentghs, which are equal to ao(00; 00)|str = 0.039 and ao(H — ; 00)|str = —0.082, respectively, the 
electromagnetic corrections are however comparable in size to the two loop strong interaction 
effects. In both cases the contributions from the counterterms fcj, estimated as mentioned above, 
are reflected by the uncertainties in these numbers and thus are seen to remain numerically 
small. Although the analysis at next-to-leading order has been carried out in the standard case, 
it is our feeling, based on the discussion at the end of Section 2, that the conclusion concerning 
the size of electromagnetic corrections to the S-wave scattering lengths will also hold within 
the generalized framework. Hopefully, an explicit calculation will confirm this expectation in 
the future. 

As a final remark, we ought to stress that the procedure followed in order to define the 
scattering length ao(H — ; 00) in Section 5 would be rather natural if we had direct experimental 
access to 7r+ — vr^ scattering. This, however, is not the case. Instead, one has to proceed along 
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indirect ways, like for instance Kg4 decays. This process, involving also leptons, has radiative 
corrections of its own, and which are only partly covered by the analysis we have presented 
here. Thus, they require a specific treatment, which should preferably be done in a systematic 
approach. This in turn means that the framework described in Sections 2 and 3 has to be 
extended, in order to accomodate the leptons. This next step, however, is beyond the scope of 
the present paper. Our purpose here was rather to assess the typical size of electromagnetic 
corrections, given a sensible, but not necessarily unique, definition of scattering lengths, and 
to exhibit, as illustrated in Figs. 2 and 3 for instance, some of the new features that emerge in 
low energy ir — ir scattering once electromagnetic interactions are taken into account. 
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Appendix 



Al. Divergent Part of the One Loop Functional 



At next-to-leading order, the generating functional is 



(A.l) 



where the integration over the fields is carried out in the one loop approximation. To evaluate 
the divergent part of the one loop functional, we transform the lagrangian of ( p.2| ) to 
Euclidian spacetime, — * C^p , and expand the fields U and around the classical solutions 
U and Afj^ of the equations of motion. 



U 



u 



A. 



U + -<M - ^< U + 

A/j, + e^. 



(A.2) 



into 



where U = u'^, and ^ is a traceless hermitian matrix. We next insert the expansion 
the action Se, and keep only the terms which are at most quadratic in the quantum fluctuations 
of the meson and of the photon fields. Collecting these fiuctuations in rj^, where A runs from 
1 to Nf + 3, ?7 = {^^, . . . , • • • , e^) we obtain, at the one loop level. 



SeU loop = lJ d'xE { r/^ + A 



a — 1 



(A.3) 



where and A are {Nj + 3) X [N] + 3) matrices. 



\xf / 



A 



(A.4) 



with 



a 



ab 



i([A„ A1[A„ A'']) + l(a{A^ A^) - \f\H,\-){H,\' 



25 



= F{([HR,A,] + ^D,Hi}jX'') , 

= -X^," = -\f{HlX^)6p^ , (A.5) 



and 



111 
if/j = u^Qru + uQlu^ , 

i^L = - . (A.6) 

The expression for a"'^ as given in Eq. (21) of was not symmetrized with respect to the 
indices a and 6, as was pointed out by the authors of Ref . p2| . Using the correctly symmetrized 



expression given above, one finds that the divergences of the counterterms -ft'15,16,17 in Eq. (35) 
of have to be replaced by 

= 3/2 + 3Z + 14^2 , S16 = -3 - 3Z/2 - , S17 = 3/2 - 3Z/2 + 5Z^ . (A.7) 

In the Feynman gauge a = 1, we thus obtain a gaussian integral for the generating 
functional, so that, up to an irrelevant constant contribution, 

-^cloneioop = ^ln(detX>) , (A. 8) 

with T)^^ = —Ti^Tif^S^^ + A^^. To renormalize the determinant we use dimensional regu- 
larization. In Minkowski spacetime, the one loop functional in c? = 4 dimensions is given by 

Zone loop = -^t^J^ J Tr V^"*" + l^^) + finite parts, (A.9) 
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where Tr means the trace in the "flavour" space t]^ and Y^^, denotes the field strength tensor 



of i;, 



Y^, = d,Y, - d,Y, + [F^, Y,]. (A.IO) 



For a 7^ 1, we obtain, in Eq. (|A.3| ), a second order differential operator of the so called 
non minimal type. Heat kernel techniques, based on pseudodifferential operator methods 



have been developped for operators of that type in the literature (see, for instance, |^ and 
references therein). They would allow, for instance, to study the gauge dependence of ^oneioopj 
but we shall not further discuss this point in the present article. 

In the Feynman gauge, the divergent part for an arbitrary number of fiavours Nf results 
from a straightforward evaluation of the trace in Eq. ( |A.9|) . If we retrict ourselves to sources 
satisfying the condition {Ql) = ( Qr ) = {Q), where Q is the constant charge matrix for Nf 
fiavours, the result becomes^ 

iTr(r^^V) + ^Tr(a2) = 

Nf Nf 
48 24 



+1 { d^'U+d'^U ) {d^U+dM ) + ( d^'U+d^U ) {d-'U+dM ) 



Nf 

-i^ ( G^d''Ud''U+ + G^^J'^U+d'^U ) 

Nf Nf 
{ G^UG^f"'U+ ) 1{ G^G^^"" + G^G^f"" ) 

Nf 1 
+^ ( d'^U^d.Uix^U + U+x) ) + ^{ d^U+d,U ){x^U + U+x ) 

o o 

N'j - A 

+ Tb^ < X^Ux^U + U^xU+x ) 



o 

( d^U+QnUd^U+QnU + d^UQLU+d^UQLU^ ) 



4 
3F2 



( dm^U+Ql + d'^U+d^UQl ) 



4 

i _ { d^^U^d^UQLU^QnU + d'^Ud.U'-QRUQLU^ ) 

+ZF^ ( d'^U+d^U ) ( QrUQlU^ )-ZF\ d^U+d^UQLU^ + d^Ud^U+QnU ){Q) 



We have written U and A instead of U and A. 
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+^ { d^U+Qnd^UQL ) + 2ZF'' { cI^U+QrU ) ( d^UQLU+ ) 

( (Ux^ + xU^)Ql + {u^x + x^u)Ql ) 
+ X^QrUQl + xQlU-'Qr ) 

+ (l + ^) { x^UQlU^QrU + xU^QrUQlU^ ) 
-ZF-" { + xU'-)Qr + ix^u + U+x)Ql ){Q) 

+ZF' { x'-U + U+x){ QrUQlU^ ) 

+^ ( d,U+[{c^^QR), Qr]U + d,U[{clQL), Ql]U^ ) 

+^ ( d^U^QnUic^QL) + d^UQLU+{di,QR) ) 

+ (2Z + 2NfZ''^ ( {QrUQlU+Y )-(2Z- 2NfZ^^ { Q\UQIU+ ) 
-SZ^F^ { QIU+QrU + QIUQlU+ ){Q) 

+ + 8Z'j F' { QrUQlU+ f-^{Ql + Ql){ QrUQlU"- ) 
+ {l + Z')F'{Ql + Qir- Z'F\ Q\-Ql)\ (A.11) 
where the covariant derivatives of the sources Ql{x) and Qr{x) are 

c'^Qi = d,Qi - t[Gl Qi], I = R,L, (A.12) 

whereas Gjj^ and Gj^^ are the field strength tensors of Gjj; and G^, respectively, 

GL = d,Gi - d^Gi - i \gL Gi] , I^R,L. (A.13) 



For Nf — 2, and after having used the equations of motion and the trace identities for 
2x2 matrices, the above expression reduces to 

^Tr(r'^^V) + iTr(a^) = 

^ ( #C/+d^C/ y + ^ {d'^U+d'^U ) { d^U+dM ) 
-^{x^U + U+Xy + li d^U+d^x + d^x-'d^U ) 

-- ( G^^UG^^^U^ )--{ G^^d^Ud''U+ + G^di^U+d^'U ) 
6 6 

+^ ( X+X ) - ^ ( G^G^'^'^ + G^^G'^^ ) + i?e(detx) 
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( d,U+d,U ){Ql + Ql) + [l-Zy'{ rU+d,U ){Qf 
+2ZF^ { d^'U+d^U ) { QrUQlU+ ) 

-^{{d^^U+QRU){d^U+QRU) + {d^^U+QLU+){d,UQLU+)) 
{ d^U+QnU ) { d^UQ^U^) - ^ {x+U + U+x) {Qr + Ql) 
-ZF' {x^U + U''x){Qf+[l + 2Z) { x^U + [/+x ) ( QrUQlU' 
+ (1-z)f^ {xU+ - Ux'-)QrUQlU+ + ix^U - U+x)QlU+QrU) 



+— { d,U+[{c^M,QRp + d,U[{clQL), Ql]U+ ) 

+^{d^U+QRU{clQL) + d^UQLU+{c^^QR) ) 

3F^ 

Y 

3 3Z 



+ (- + 3Z + 12^2) F' { QrUQlU+ QrUQlU+ ){QI + QI 

- (3Z + 12^2) F' { QrUQlU^ ){Qf+(^^-^ + z'y'{Ql + Q 



2 \2 
L I 



CiZ 



^\^-:--2Z^]F"{Q\^Ql){QY- Z^F' {Ql^Qlf^ AZ^F' ( Q )^ (A. 14) 



A2. Loop Functions 



We recall the definition of the function Jpq{s) = Jpq{s) — Jpq(O), with 



JpQis) = / 7^ 



d^ 1 

(2^0^ -Ml> + ie)[iq - - + ie] ' 



and s —p^. For s > {Mp + Mq)'^ and for d — A, one finds 



327rVpQ(s) = 2 + 



In 



s 

A|>o(s) 



Ml ApQ Ml 



s - \lQ{s)f - A 



PQ 



+ 2i7r-^^ 



with 



and 



^PQ = M| + , ApQ = Ml -Ml, 
Xpq{s) = \{s,Ml,Ml) . 



(A.15) 



(A.16) 

(A.17) 
(A.18) 
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For the threshold expansions of the amphtudes or cross sections, the following expressions have 
been used, 



i?eJoo(4M2±) 
i?e J+_(4M^o) 
i?e J+o(-A,) 



1 



2 A2 



1 + 



Ml, 3 M^o 
A. 1 A! 



M^o 3 M\ 



1 A, 



i?e J+o((M^± +M^o)2 



967r2 M^o 
1 



^_ 3 A, 



5 M^, 
1 A2 



48 M\ 



1 A2 



3847r2 M^o 



(A.19) 



where only the contributions up to order 0(A^/M^o) have been kept. 

As far as the function defined in Eq. (|4.12| ) is concerned, introducing two Feyn- 

man parameters and performing standard manipulations leads to the following integral repre- 
sentation 



^+-7('^) 



327r2 Jo 



dx 



dy 



d 



f{x) dy 



In 



y'^fi^) + (1 - y)"^: 



+ Oim, 



7/ ? 



where 



f{x) = Ml± - x{l - x)s -ie 



(A.20) 



(A.21) 



For s < 0, the roots of f{x) are real and lie outside of the interval [0, 1 ], so that the integration 
is straightforward. 



1 



4Li2 



1 — cr\ TT 



1 + (T 



a-1 



In 



In' 



+ 21n(cr) 



In 



a - 1 
a+1 



(A.22) 



The dilogarithm or Spence function is defined as usual, 

Li2(a;) 



1 1-t 



dt . 



(A.23) 



The expression ([4.13 ) for s > 4:Ml± follows by analytic continuation with the ie prescription. 
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A3. Radiative vr — vr Scattering 



Let 



s - 28^ - AMI, 



S - 2Sy^ ■ 



(A.24) 



Then the function F{s;AE), resulting from the integration over the energy of the emitted 
photon, is given by the indefinite integral 



2 \2 ■ 



/ISl^ Hi- 
— [{s-28^s-Ml,ff{£,s)-{s 

- i(s - 2AEv^)(6M2„ + 2AEv^ - 3s)f{AE, s) 

~{s - Ml,f In I [ s - 2AE^s - 4M^o ] [ 1 + /(AE, s) ] + 4M^oAE [ 1 + 2/(A£;, s) ] | 

(A.25) 



s - AMI, 
In Us - AM\ 



s' - AsMfo + SMlo] 



{s-2AEy^)[l 



'—^f{AE,s)]-2Ml, 



+ Cst . 



For the threshold expansion of the total cross section, we need 

F{AMl^-^) - F{AMl^-Q) = 

M^± 

l2Ml±{Mlo - 2Ml±) + {AMl± - M^of In 



'AM^ 
. Ml, 



^i:^ ilQMt^ _ iQMl^Ml, + 2,Mt.) In { '^^^ ^M^^V^^ 



M^ 

18(ln2-l)M^o +■ 



(A.26) 



Since the difference F{AM^±; A^/ M^±) - F{AM^±;0) already appears multiplied by e^A^ in 
the expression ( ^.141 ) of a^^'^^'' at threshold, we have dropped, in the last line of (A.26), all 
the contributions that vanish for Att 0. 
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